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ABSTRACT 

We discuss the interpretation of transient, unpulsed radio emission detected from the 
unique pulsar/Be-star binary system PSR B1259— 63. Extensive monitoring of the 1994 
and 1997 periastron passages has shown that the source flares over a 100-day interval 
around periastron, varying on time-scales as short as a day and peaking at 60 mJy (~ 100 
times the apastron flux density) at 1.4 GHz. Interpreting the emission as synchrotron 
radiation, we show that (i) the observed variations in flux density are too large to be 
caused by the shock interaction between the pulsar wind and an isotropic, radiatively 
driven, Be-star wind, and (ii) the radio emitting electrons do not originate from the 
pulsar wind. We argue instead that the radio electrons originate from the circumstellar 
disk of the Be star and are accelerated at two epochs, one before and one after periastron, 
when the pulsar passes through the disk. A simple model incorporating two epochs of 
impulsive acceleration followed by synchrotron cooling reproduces the essential features 
of the radio light curve and spectrum and is consistent with the system geometry inferred 
from pulsed radio data. 



1. Introduction 

The binary pulsar PSR B1259— 63 is in a highly 
eccentric (e ~ 0.87), 3.5-yr orbit around a 10th mag- 
nitude B2e star, SS 2883, of radius i?* ~ 6 R© (John- 
ston et al. 1992, 1994). The system is unique because 
PSR B1259— 63 is the only radio pulsar orbiting a 
main-sequence companion with a circumstellar disk. 
Observations of the Ha emission line show that op- 
tically emitting material in the Be-star disk extends 
to at least 20 i?* (Johnston et al. 1994), and pulsar 
timing data suggest that the disk is steeply inclined 
to the plane of the binary orbit (Wex et al. 1998). 
The pulsar's orbital radius at periastron is 23 i?*. 
Melatos, Johnston & Melrose (1995) proposed a wind- 
disk model of the system to explain orbital variations 
of the flux density, polarisation, rotation measure and 
dispersion measure of the pulsed emission near pe- 
riastron (Johnston et al. 1996). In the model, the 
eclipse of the pulsed emission between T — 20 and 
T + 16 (where T denotes the epoch of periastron) is 
attributed to free-free absorption in the disk. 

X-ray observations of PSR B1259— 63 reveal sig- 
nificant (unpulsed) emission throughout the orbit, 
with the flux density at periastron ~ 10 times higher 
than at apastron (Hirayama et al. 1996). Modeling 
suggests that the X-rays are synchrotron radiation 
from rclativistic pulsar wind leptons accelerated at 
the shock between the pulsar wind and Be-star out- 
flow (Tavani & Arons 1997). 

We discuss the origin of the transient, unpulsed 
radio emission detected from PSR B1259— 63 from 
T-22 until T +100 (Johnston et al. 1999). The radio 
emission, like the X-rays, is likely to be synchrotron 
radiation produced by the shock interaction between 
the pulsar wind and Be-star outflow. However, we ar- 
gue that the interaction between the pulsar wind and 
an isotropic, radiatively driven, Be-star wind cannot 
explain the magnitude and rapidity of the variations 
in the radio light curve. Furthermore, we show that 
the radio emitting electrons, unlike the X-ray elec- 
trons, do not originate from the pulsar wind. We 
therefore suggest that the unpulsed radio emission re- 
sults from the acceleration of electrons in the Be-star 
disk. The pulsar-disk encounter is likely to result in a 
strongly time-dependent, magnetized, hydrodynami- 
cal interaction forming a cometary pulsar wind neb- 
ula surrounded by a shock front where acceleration 
of particles is followed by a combination of adiabatic, 
synchrotron and inverse Compton cooling. However, 



the data so far obtained from this system are insuffi- 
cient to constrain a model of the physics at this level 
of detail. Here we propose a simple model for the ra- 
dio light curves which reproduces the main features 
of the data and identifies the essential physics. 

2. Observations 

Two periastron passages of PSR B1259— 63 have 
been observed in detail, the first on 1994 Jan 9 and the 
second on 1997 May 29. Unpulsed radio emission was 
detected at frequencies between 0.84 and 8.4 GHz. 
as described by Johnston et al. (1999). Figure 4 of 
that paper shows the radio light curve at five frequen- 
cies during the 1997 periastron passage (the observed 
1.4 GHz light curve is reproduced here in Figure lb). 
Emission appeared suddenly at T — 22; prior to this 
the 3-cr upper limit on the unpulsed emission is < 0.5 
mJy at frequencies between 1.4 and 8.4 GHz. For the 
next few days, rapid variations were superimposed on 
an overall increase in the flux density S at all fre- 
quencies, peaking at S ~ 30 mJy at 1.4 GHz near 
T — 8. There followed a relatively slow decline before 
the flux density again began increasing around T + 8, 
peaking a second time near T + 18, with S ~ 60 mJy 
at 1.4 GHz. The flux density then declined steeply 
until around T + 30, when the rate of decay slowed to 
roughly that observed after the first peak near T — 8. 
This slower decline continued until the last observa- 
tions on T+47. The 1994 data are similar and showed 
detectable emission until T + 100. 

At most times the radio frequency spectrum is well 
fitted by a power law S oc v~ a of spectral index 
a 0.5, consistent with synchrotron radiation from 
relativistic electrons. Importantly, there is evidence 
in both the 1994 and 1997 data that the spectrum 
steepens after periastron (Figure 5 of Johnston et al. 
1999). There is no compelling evidence for free- free 
absorption at any epoch between T — 22 and T + 47 
(Johnston et al. 1999), even at periastron when the 
pulsed emission is totally eclipsed by the Be-star disk. 

3. Origin of the radiating electrons 

The relativistic pulsar wind is enclosed within a 
bubble as it ploughs through the material around the 
Be star. Momentum flux balance between the pul- 
sar wind and Be-star outflow determines the size and 
shape of the bubble, and particles are shock acceler- 
ated at the boundary between the two flows (Melatos, 
Johnston & Melrose 1995). 
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The Be-star outflow consists of a radiation-driven 
polar wind and a more dense equatorial disk re- 
sponsible for the bulk of the mass loss (Dougherty 
1994). Both components are often characterised as 
simple power laws in radial distance r s , with num- 
ber density and speed given by n s = n (i?*/r s )£ and 
Vs = vo(r s / R*)^~ 2 , where £ = 2 corresponds to a 
constant speed and £ = 4 corresponds to constant 
ram pressure. If the outflow is highly conducting so 
that the magnetic field is frozen into the plasma, then 
simple models imply B s oc r^~^ for 2 < £ < 3. The 
polar wind is generally thought to reach a terminal 
speed v w <~ 1000 kms -1 within a few stellar radii 
(e.g. Friend & McGregor 1984); we therefore approx- 
imate it as a constant speed outflow with £ w = 2. 
In contrast, the disk density falls off rapidly with ra- 
dius, as rotating material gradually accelerates out- 
ward; one typically has £d~4, v d0 ~ 5 kms -1 and 
a half opening angle of a d ~ 5 — 15° (Waters, Cote 
& Lamers 1987). The disk and wind number densi- 
ties are related by n d0 /n w0 = (v w /v d0 )(M d / M w )(l - 
sin a d ) / sin a d . For a d = 10° this implies a density 
ratio of n d0 /n wn <~ 10 3 (M d /M w ). Typical values 
of M ~ 2 x 10~ 8 Meyr 1 and M d /M w ~ 10 yield 
n d i)~4x 10 18 m~ 3 . 

3.1. Flux scaling for polar wind interaction 

If the spin-down luminosity of the pulsar, E p , goes 
entirely into the relativistic wind then the ram pres- 
sure at a distance r p from the pulsar is E p /{^ircr p ). 
The ram pressure of the Be-star polar wind is 
m p n s (r s )v s (r s ) . The contact discontinuity between 
the two flows has its apex at the point where the 
ram pressures balance, provided the magnetic pres- 
sures can be ignored. This point is denoted by 
r p = D — r s = L where D is the pulsar-Be star sepa- 
ration, L/D = ^/rj/{l + y/rj) and r\ = E p /(M w cv w ) is 
the ratio of the momentum flux in the pulsar wind to 
that in the Be-star polar wind. The ratio L/D is con- 
stant over the pulsar orbit, although D and L both 
vary. 

Now suppose that the unpulsed radio flux is asso- 
ciated with the shock interaction between the pulsar 
wind and Be-star polar wind. Synchrotron theory im- 
plies that the flux density scales as 

S <xnl 3 B a+1 v- a (1) 

where n is the number density of radiating electrons, 
I is the source size and B is the magnetic field. As- 
suming I oc L and that the other source parameters 



scale as their values at the contact surface apex, one 
obtains I oc D, n oc D~ 2 and B oc D^ 1 , implying 
S oc D~ x / 2 for a = 0.5 (irrespective of whether the ra- 
diating electrons originate from the pulsar wind or Be- 
star polar wind). Since D varies by a factor of <~ 15 
around the orbit, this argument implies that S should 
be ~ 4 times higher at periastron than at apastron — 
consistent with the sparse X-ray data, but far smaller 
than the 1.4 GHz radio variation from 5<0.5 mJy 
at apastron to S ~ 60 mJy near periastron. Fur- 
thermore, in the interval near periastron, the above 
scaling yields S(T - 27) - S(T)/y/3 ~ 20 mJy, yet 
at T — 27 there is no detectable unpulsed radio emis- 
sion at any frequency. We conclude that the emission 
cannot be attributed solely to the interaction between 
the pulsar wind and Be-star polar wind. 

3.2. Disk or pulsar electrons? — 
Lorentz factor 

The pulsar orbit is eccentric and steeply inclined 
with respect to the Be-star disk, so the pulsar wind 
interacts strongly with the disk near periastron, po- 
tentially generating unpulsed radio emission. We now 
investigate whether the radiating electrons in this sce- 
nario originate from the pulsar wind or the disk. 

The magnetic field in the emission region is the 
sum of the contribution from the pulsar, B p (L) = 
1.5 x 10~ 3 i?*/L tesla (estimated from P, P and the 
assumption B p oc r p 1 in the pulsar wind), and the 
contribution from the Be star. The field in the Be-star 
polar wind satisfies B S (D-L) < 10~ 2 i?„/(Z) - L) tesla 
(see Barker 1987 for an upper limit on surface fields 
of Be stars), while estimates of the field in the disk 
are so uncertain as to be of little value. For L<D/2, 
one finds B ~ 10~ 4 T at D ~ 50i?» when the ra- 
dio peaks occur, and hence B p (L) ~ B S (D — L). 
This value of B is consistent with that obtained inde- 
pendently from synchrotron equipartition arguments, 
which yield B - 10~ 4 T for S(l GHz) - 50 mJy 
and a = 0.5, assuming a source distance of 1.5 kpc, a 
spherical source of radius 25 

The Lorentz factor required to produce synchrotron 
emission at 1 GHz in such a field is 7 ~ 50(r p /50i?») 1 / 2 
- five orders of magnitude lower than the likely 
Lorentz factor 7 ~ 10 6 of unaccelerated leptons in the 
pulsar wind (Kennel & Coroniti 1984; Arons & Gal- 
lant 1994). We conclude that the pulsar- wind elec- 
trons arc too energetic to be responsible for the un- 
pulsed radio emission observed. On the other hand, 
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non-rclativistic electrons from the Be-star disk and/or 
polar wind could certainly be accelerated to the rel- 
atively modest energy required at shock(s) produced 
by the interaction between the pulsar wind and Be- 
star disk (e.g. Kirk 1994) as they are in binaries com- 
prising early-type stars (Eichler & Usov 1993). 

3.3. Disk or pulsar electrons? — 
Number density 

The number density of radiating electrons can be 
estimated straightforwardly for an assumed source 
size and magnetic field, since from synchrotron the- 
ory n oc S l~ 3 B~ 3 / 2 . The source size I is at least 
the standoff distance L, otherwise synchrotron self- 
absorption would cause a detectable low-frequency 
turnover below 1 GHz (Johnston ct al. 1999). If the 
emission results from the interaction between the pul- 
sar wind and Be-star disk, then the source is likely to 
be some kind of sheath around the contact disconti- 
nuity, extending downstream in a cometary tail. A 
spherical sheath with thickness x L, where £ is a nu- 
merical factor less than 1, yields n ~ x^ 1 10 10 m~ 3 
for the parameters used above. 

The number density of the pulsar wind can be es- 
timated straightforwardly from the spindown lumi- 
nosity E p = 8.3 x 10 28 W and wind Lorentz factor 
7 <~ 10 6 , yielding n p (r p ) = E v / A-K"fm e c 3 r p and hence 
n p (L) - 2 x 10 4 m~ 3 for L ~ 25.R*. This is at least 5 
orders of magnitude smaller than the estimated num- 
ber density required to produce the observed radio 
flux density. While the estimates clearly allow consid- 
erable latitude, and a shock may compress the wind 
by a factor of order 10, it is very difficult to argue 
away this discrepancy. We conclude that the pulsar- 
wind electrons are too few in number to be respon- 
sible for the unpulscd radio emission observed. On 
the other hand, even if the number density of the Be- 
star disk falls off as steeply as n<j oc r~ 4 , one finds 
nd(25i?*) <~ 10 13 m~ 3 , and acceleration of just 1CP 3 
of the electrons in the disk can account for the ob- 
served emission. 

4. Model of the radio light curve 
and spectrum 

Johnston ct al. (1999) and Wex et al. (1998) have 
argued that the Be-star disk is inclined with respect to 
the orbital plane in such a way that the pulsar passes 
through the disk twice in each orbit, once before peri- 
astron (crossing the disk around T— 18 at Z?i ~ 40i?*) 



and once afterwards (T+ 13 at D 2 ~ 33i?»). We pos- 
tulate that the radiating electrons are created dur- 
ing time intervals T\ and T 2 at the two disk cross- 
ings, when the shock interaction between the pulsar 
wind and the disk (and hence electron acceleration) 
is at its peak, and nowhere else. Ti and T 2 are some- 
what greater than the orbital crossing time of the disk 
(~ 3 days with = 5° according to Johnston ct al. 
1999) because the pulsar wind interacts appreciably 
with the disk surface for a few days before the pulsar 
enters and after it leaves. 

The two resulting electron populations are as- 
sumed to evolve with time purely under the ac- 
tion of synchrotron losses. The loss times at the 
two crossings are given in terms of constants 61,2 = 
C(B\ 2 /2/j, ), where is the local magnetic field, 
C = 327rrg/(9TOgC 3 ), and r e is the classical elec- 
tron radius. The number densities per unit energy 
^i,2{e, t) then obey 

g^jM = hl2 ^ Nla[£it )] + Q h2 ( £ ,t), (2) 

where Qi t 2(e,t), the injection rate at each disk cross- 
ing, is a power law in energy during the times T\ and 
T 2 and zero otherwise. Note that the evolution of 
the accelerated electrons depends only on conditions 
at the acceleration site; it is unrelated to the subse- 
quent orbital position of the pulsar. Emission from 
the first disk crossing persists as the pulsar rounds 
the Be star (behind the disk) and is still present when 
the pulsar re-enters the disk after periastron, so the 
observed emission is the sum of that from the two 
electron populations. 

The lack of evidence of a low-frequency spectral 
turnover implies that there is little absorption of the 
unpulsed radio emission (Johnston et al. 1999), even 
when the pulsar is at periastron and is obscured be- 
hind the dense Be-star disk. This strongly suggests 
that the radio-emitting electrons are left behind by 
the pulsar in a ram-pressure-confined bubble at the 
point where the pulsar crosses the disk, and that 
therefore their radio emission is not obscured by the 
bulk of the disk (Ball et al. 1998). We therefore ig- 
nore propagation effects on the radio emission. 

The proposed model does not include the effects 
of adiabatic expansion in the evolution of Ni^{e,t). 
Melatos, Johnston & Melrose (1995) showed that 
when the pulsar is deep within the Be-star disk, the 
momentum flux in the pulsar wind is insufficient to 
disrupt the disk and so the pulsar wind is contained 
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within a bubble. It follows that in a model which as- 
sumes that the bubble of radiating electrons is spheri- 
cal, the source does not expand and so adiabatic losses 
do not occur. Electrons accelerated during each disk 
crossing are advected away from the shock apex in 
a cometary backflow (Melatos, Johnston & Melrose 
1995), but the observed rate of decline of the unpulsed 
radio emission indicates that the resulting adiabatic 
losses must occur on a timescale of a week or longer. 
In this Letter we show that a model involving only 
injection during the pulsar disk crossings and syn- 
chrotron losses is in qualitative agreement with the 
observations. 

Inverse Compton losses may be important in this 
system (Tavani & Arons 1997; Kirk, Ball & Skjaeraasen 
1999). At the energies of the radio-emitting electrons 
these losses have the same form as synchrotron losses, 
and so the evolution of the electron spectrum will 
still be described by an equation of the form of (g). 
Outside of the disk the energy density of target pho- 
tons from the Be star is comparable to the energy 
density in the magnetic field at the shock between 
the pulsar wind and the Be-star wind (Kirk, Ball & 
Skjaeraasen 1999), in which case the inverse Comp- 
ton and synchrotron loss rates would be compara- 
ble. However, the density of target photons available 
to radio-emitting electrons originating in the disk is 
likely to be low because of screening by the disk (Ta- 
vani & Arons 1997), and so inverse Compton losses 
are likely to be less important than synchrotron losses 
in the model for the radio emission proposed here. 
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Fig. 1. — (a) Model flux density F{y) in arbitrary 
units as a function of time (solid line) , and at a fre- 
quency v' — 1.44 v (dashed line), (b) Observed flux 
density from PSR B1259-63 at 1.4 GHz (from Figure 
4 of Johnston et al. 1999). 

Assuming that the synchrotron emissivity can be 
approximated by a delta-function in frequency, the 



observed emission at a given frequency v originates 
from electrons with characteristic energies e± and £ 2 
at the two acceleration sites. Since v cx Bi 2 s\ 2 it fol- 
lows that in general E\ 7^ £2- The total number den- 
sity per unit energy of electrons emitting characteris- 
tically at v is then Ni(ei,t) + N 2 {s2, t). When adia- 
batic losses (and therefore expansion) are neglected it 
follows from (|lj) that the that the emitted flux density 
is proportional to the number density of electrons ra- 
diating at v. The solid line in Figure la is the model 
flux density F(v) obtained in this way, plotted in ar- 
bitrary units as a function of time. The curve is ob- 
tained by solving equation (||) for T\ = T 2 = 10 d, 
l/(6i£i) = 72 d, l/(&2£2) = 9d, and source functions 
of the form Qi i2 (e,t) cx (£/£ )~ a {e min < e < e max , 
a = 2.2). The dashed curve shows the model flux 
density for a frequency v 1 = l.AAv. Since the pul- 
sar is closer to the Be star during the post-periastron 
disk crossing than during the pre-periastron cross- 
ing, ram pressure balance implies L 2 < L\; however, 
the choice L1/L2 = 4 (= B2/B1) in Figure la is 
purely illustrative because we have no precise knowl- 
edge of the geometry of the pulsar-disk interaction. 
The synchrotron-loss times l/(£>i,2£i,2) are roughly 
the values expected for B ~ 10 -4 T and v ~ 1 GHz; 
the injection intervals T\ 2 are roughly twice the pul- 
sar disk crossing times; the relative magnitudes of 
Qi and Q2 are chosen arbitrarily to be similar. The 
curves are approximately piecewise-linear, a charac- 
teristic of solutions to equation ^ for a w 2. The 
proposed model reproduces many of the essential fea- 
tures of the light curves in Figure 4 of Johnston et al. 
(1999), from which the 1.4 GHz light curve is repro- 
duced in Figure lb. In particular, the flux density is 
higher after T + 12 because it is the sum of the con- 
tributions from two accelerated electron populations. 

The following features of the light curves in Fig- 
ures la and lb should be noted, (i) The first rise is 
associated with the pre-periastron disk crossing and 
lasts for a time T\ , because the synchrotron- loss time 
l/(fri£i) is much greater than T\. (ii) The gradual 
decay after T — 8, which continues until T + 54 (solid 
curve), occurs because the number of electrons per 
unit time cooling down to energy e\ from e > E\ 
is slightly smaller than the number with energy e\ 
cooling to lower energies. (Exact balance occurs for 
a = 2.) (iii) The second rise is associated with 
the post-periastron disk crossing and lasts for the 
synchrotron-loss time l/(&2£2), because T 2 is larger 
than l/(&2£2)- (iv) The plateau at the principal max- 
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imum lasts for a time Ti — l/(&2£2) — Id, which 
is brief because the injection and loss times for the 
post-periastron disk crossing have been chosen to be 
almost equal, [cf. (ii)] (v) Following this brief plateau, 
7V 2 (e2,i) decreases to zero over a time 1/(62^2) equal 
to the rise time in (hi). Since ^V 2 (e 2 ,t) decreases to 
zero relatively rapidly, its contribution drops below 
that of the more slowly decaying N\(ei,t) and the 
subsequent evolution is therefore an extension of (ii) . 
Note that the pre-periastron rise time (i) is deter- 
mined by injection (which lasts longer than the loss 
time) , whereas the post-periastron rise time (iii) is de- 
termined by the loss time (which is shorter than the 
second injection episode). These rise times happen 
to be similar because we chose parameters such that 
Ti ~ 1/(6262). 

The changing separation of the solid and dashed 
curves in Figure la indicates that the effective spec- 
tral index a of the composite electron distribution is 
not constant. By crudely relating a to a (S oc v~ a ) 
via a = (a — l)/2, we predict that the spectrum: (i) 
has a ~ 0.7 and does not change significantly from 
switch-on at T — 18 until the short plateau around 
T + 20; (ii) then steepens markedly with a increasing 
to ~ 1.1 before flattening again to <~ 0.7 at T + 30; 
(iii) steepens smoothly thereafter as the emission at 
higher frequencies decays more rapidly. The model 
behavior is broadly in accord with the data, which 
reveal a systematic steepening of the spectrum with 
the spectral index changing by ~ 0.4 between T + 15 
and T+30 (Figure 5 of Johnston et al. 1999), followed 
by a less pronounced flattening. 

The model makes a clear, testable prediction. Wc 
see from Figure la that the slow decay after the first 
disk crossing (which continues well past the second 
disk crossing to the end of the observation period 
T+47) ends sharply; the model flux density decreases 
to zero over a time comparable to the initial rise, 
T\ = 10 d, instead of maintaining its very gradual 
downward slope all the way to zero. Although obser- 
vations extending to T+ 100 were made at 0.84 GHz 
following the 1994 periastron, these single-frequency 
data are inconclusive as to the existence of such a 
feature. It will be possible to test this prediction dur- 
ing the next periastron passage; the frequency depen- 
dence of the onset of this rapid decline should also be 
observable. 

The model makes no attempt to explain the short- 
est time scale variations in the radio emission. Partic- 
ularly puzzling is the flat spectrum 'blip' which was 



observed at T — 21 in 1997. We speculate that this 
feature is the result of the pulsar nebula splashing 
into the Be-star disk. Such effects are more likely to 
have observable consequences the first time the pulsar 
enters the disk because the system geometry is such 
that we have an unobscured view of the pre-periastron 
impact (see Figure 6 of Johnston et al. 1999), but not 
of the post-periastron impact. 

5. Conclusions 

We have argued that the transient, unpulsed radio 
emission from the binary system PSR B1259— 63 is 
synchrotron emission associated with the shock inter- 
action of the relativistic pulsar wind and the dense cir- 
cumstellar disk of the Be star. The radiating electrons 
originate from the disk, not a spherically-symmetric 
pulsar wind with the parameters required to model 
the X-ray data. 

We have presented a simple physical model for the 
evolution of the unpulsed radio emission, based on 
impulsive acceleration of electrons followed by syn- 
chrotron losses, which reproduces the salient features 
of the radio light curves and spectra observed during 
the 1994 and 1997 periastron passages. The model 
also makes a testable prediction of a rapid frequency- 
dependent decay in the radio flux several weeks after 
the principal flux density peak, to be searched for 
during observations of the next periastron passage in 
October 2000. 
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